Spin-Orbit Coupling, Antilocalization, and Parallel Magnetic Fields in Quantum Dots 
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We investigate antilocalization due to spin-orbit coupling in ballistic GaAs quantum dots. Antilo- 
calization that is prominent in large dots is suppressed in small dots, as anticipated theoretically. 
Parallel magnetic fields suppress both antilocalization and also, at larger fields, weak localization, 
consistent with random matrix theory results once orbital coupling of the parallel field is included. 
In situ control of spin-orbit coupling in dots is demonstrated as a gate-controlled crossover from 
weak localization to antilocalization. 



The combination of quantum coherence and electron 
spin rotation in mesoscopic systems produces a number 
of interesting and novel transport properties. Numerous 
proposals for potentially revolutionary electronic devices 
that use spin-orbit (SO) coupling have appeared in recent 
years, including gate-controlled spin rotators jlj as well 
as sources and detectors of spin-polarized currents Q . It 
has been predicted that the effects of some types of SO 
coupling will be strongly suppressed in small OD systems, 
i.e., quantum dots || ||, This suppression as well 
as overall control of SO coupling will be important if 
quantum dots are used to store electron spin states as 
part of a future information processing scheme. 

In this Letter, we investigate SO effects in ballistic- 
chaotic GaAs/AlGaAs quantum dots. We identify the 
signature of SO coupling in ballistic quantum dots to be 
antilocalization (AL), leading to characteristic magneto- 
conductance curves, analogous to known cases of disor- 
dered ID and 2D systems § % § § ^ Q. AL is 
found to be prominent in large dots and suppressed in 
smaller dots, as anticipated theoretically J|, |J, ■ Results 
are generally in excellent agreement with a new random 
matrix theory (RMT) that includes SO and Zeeman cou- 
pling H . Moderate magnetic fields applied in the plane of 
the 2D electron gas (2DEG) in which the dots are formed 
cause a crossover from AL to weak localization (WL). 
This can be understood as a result of Zeeman splitting, 
consistent with RMT || . At larger parallel fields WL is 
also suppressed, which is not expected within RMT. The 
suppression of WL is explained quantitatively by orbital 
coupling of the parallel field, which breaks time-reversal 
symmetry [jl2|. Finally, we demonstrate in situ electro- 
static control of the SO coupling strength by tuning from 
AL to WL in a dot with a center gate. 

It is well known that in mesoscopic samples coherent 
backscattering of time-reversed electron trajectories leads 
to a conductance minimum (WL) at B = in the spin- 
invariant case, and a conductance maximum (AL) in the 
case of strong SO coupling jfjj. In semiconductor hct- 
erostructures, SO coupling results mainly from electric 
fields |l| (appearing as magnetic fields in the electron 
frame) leading to momentum dependent spin precessions 
due to crystal inversion asymmetry (Dresselhaus term 
p^|) and heterointerface asymmetry (Rashba term [fl5||). 



SO coupling effects have been pre viously measured us- 
ing AL in GaAs 2DEGs §, |, [Uj and other 2D het- 
erostructures (yj. Other means of measuring SO cou- 
pling in heterostructures, such as from Shubnikov-de 
Haas oscillations |l6| and Raman scattering spectroscopy 
|l7f are also quite developed. SO effects have also been 
reported in mesoscopic systems (comparable in size to the 
phase coherence length) such as Aharonov-Bohm rings, 
wires, and carbon nanotubes jL8j. Recently, parallel field 
effects of SO coupling in quantum dots were measured 
p9| , p0| . In particular, an observed reduction of conduc- 
tance fluctuations in a parallel field [ p0[ was explained by 
including SO effects [Q, ||, leading to an important ex- 
tension of random matrix theory (RMT) to include new 
symmetry classes associated with SO and Zeeman cou- 
pling g. 

This RMT addresses quantum dots coupled to two 
reservoirs via N total conducting channels, with N 1. 
It assumes (7,ez) <C Et, where 7 = NA/(2tt) is the 
level broadening due to escape, A is the mean level spac- 
ing, ez = 9[J-bB is the Zeeman energy and Et is the 
Thouless energy (Table I). Decoherence is included as 
a fictitious voltage probe || with dimensionless de- 
phasing rate N v — h/(AT v ), where t v is the phase co- 
herence time. SO lengths A 1,2 along respective principal 
axes [110] and [110] are assumed (within the RMT) to be 
large compared to the dot dimensions along these 
axes. We define the mean SO length \ so — \f\X~i\~2\ 
and SO anisotropy v so = \J\\i/\2\- SO coupling in- 
troduces two energy scales: e s f = k±Et(LiL 2 / X 2 ) 2 , 
which represents a spin-dependent Aharonov-Bohm-likc 
effect, and ejj° ~ ((^i/Ai) 2 + (L 2 / X 2 ) 2 )e s f , providing 
spin flips. AL appears in the regime of strong SO cou- 
pling, (e™, ejj°) 3> 7, where 7 is the total level broadening 
7 = (i + h/Ty). Note that large dots reach the strong SO 
regime more readily (i.e., for weaker SO coupling) than 
small dots. Parameters A so , t v , and k± (a dimensionless 
parameter characterizing trajectory areas within the dot) 
are extracted from fits to dot conductance as a function 
of perpendicular field, B±. The asymmetry parameter, 
v so , is estimated from the dependence of magnetocon- 
ductance on parallel field, By . 

The quantum dots are formed by lateral Cr-Au de- 
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pletion gates defined by electron-beam lithography on 
the surface of a GaAs/AlGaAs heterostructure grown in 
the [001] direction. The 2DEG interface is 349 A be- 
low the wafer surface, comprising a 50 A GaAs cap layer 
and a 299 A AlGaAs layer with two Si <5-doping layers 
143 A and 161 A from the 2DEG. An electron density of 
n — 5.8 x 10 15 m~ 2 [|2| and bulk mobility fi ~ 24m 2 /Vs 
(cooled in the dark) gives a transport mean free path 
4 - 3 fim. This 2DEG is known to show AL in 2D 0. 
Measurements were made in a 3 He cryostat at 0.3 K us- 
ing current bias of 1 nA at 338 Hz. Shape-distorting gates 
were used to obtain ensembles of statistically indepen- 
dent conductance measurements [^3| while the point con- 
tacts were actively held at one fully transmitting mode 
each (N = 2). 

Figure 1 shows average conductance (g) , and variance 
of conductance fluctuations, var(g), as a function of B± 
for the three measured dots: a large dot (^4 ~ 8/xm 2 ), a 
variable size dot with an internal gate (A ~ 5.8 /im 2 or 
8 /im 2 , depending on center gate voltage), and a smaller 
dot (1.2 /j,m 2 ). Each data point represents ~ 200 inde- 
pendent device shapes. The large dot shows AL while the 
small and gated dots show WL. Estimates for A so , t v and 
k±, from RMT fits are listed for each device below the 
micrographs in Fig. 1 (see Table I for corresponding ej_ 
and eii). When AL is present (i.e., for the large dot), es- 
timates for A so have small uncertainties (±5%) and give 
upper and lower bounds; when AL is absent (i.e., for the 
small and gated dots) only a lower bound for A so (—5%) 
can be extracted from fits. The value A so ~ 4.4 /ira is 
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FIG. 1: Average conductance (g) (squares) and variance of con- 
ductance var(g) (triangles) calculated from ~ 200 statistically inde- 
pendent samples (see text) as a function of perpendicular magnetic 
field B± for (a) 8.0 /tm 2 dot (b) 5.8 /tm 2 center-gated dot and (c) 
1.2 /im 2 dot at T = 0.3 K, along with fits to RMT (solid curves). 
In (b), the center gate is fully depleted. Vertical lines indicate the 
fitting range, error bars of (g) are about the size of the squares. 
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FIG. 2: (a) Difference of average conductance from its value at 
large B±, Sg(B±, -By), as a function of B± for several By for the 
8.0 /im 2 dot at T = 0.3 K (squares) with RMT fits (curves), (b) 
Sensitivity of <5g(0, By) to v ao for the 8.0 /mi 2 dot, 1 < v ao < 2 
(shaded), v so = 1.4 (solid line) and u so = 0.8 (dashed line) (c) 
Sg(0, By ) (markers) with RMT predictions (dashed curves) and one 
parameter (solid curves) or two parameter fits (dotted curves) using 
RMT including a suppression factor due to orbital coupling of By , 
see text. 

consistent with all dots and in good agreement with AL 
measurements made on an unpatterned 2DEG sample 
from the same wafer JhJ . 

Comparing Figs. 1(a) and 1(c), and recalling that all 
dots are fabricated on the same wafer, one sees that AL is 
suppressed in smaller dots, even though A so is sufficient 
to produce AL in the larger dot. We note that these dots 
do not strongly satisfy the inequalities L/X so < 1,JV> 
1, having N — 2 and L/X so = 0.64 (0.34) for the large 
(small) dot. Nevertheless, Fig. 1 shows the very good 



A 


A 


Td 


E T /A 


eI°/A 


ejj°/A 


Oi, a 2 


b 2 


/tm 2 


fieV 


ns 








(ns)"^- 2 


(ns)-^" 6 


1.2 


6.0 


0.35 


33 


0.15 


0.04 


6.6, 6.6 


0.24 


5.8 


1.2 


1.7 


73 


0.32 


0.33 


3.2, 


140 


8 


0.9 


2.3 


86 


3.6 


3.1 


1.4, 0.9 


3.7 



TABLE I: Dot area A = L1L2 (130 nm edge depletion); spin- 
degenerate mean level spacing A = 27r7i 2 /m*A (m* = 0.067m e ); 
dwell time r d = h/(NA); Thouless energy E T = hvp/y^A; e s f/A 
and ejj°/A for the fits in Fig. 1; B 2 coefficients cti and ct2 from one 

and two parameter fits; _B 6 coefficient 62 from two parameter fit, 
see text. 
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FIG. 3: (a) Difference of average conductance from its value at 
large £?x, Sg(B±,0), for various temperatures with Bn = for the 
8.0 /im 2 dot (squares), along with RMT fits (solid curves), (b) Spin- 
orbit lengths Xso (circles) and phase coherence times t v (triangles) 
as a function of temperature, from data in (a). 

agreement between experiment and the new RMT. 

We next consider the influence of a parallel magnetic 
field on average magnetoconductance. In order to apply 
tesla-scale Bn while maintaining subgauss control of B±, 
we mount the sample with the 2DEG aligned to the axis 
of the primary solenoid (accurate to ~ 1°) and use an 
independent split-coil magnet attached to the cryostat 
to provide B± as well as to compensate for sample mis- 
alignment p(J . Figure 2 shows plots of the deviation of 
the shape-averaged conductance from its value at B± 3> 
4>o/A (i.e., with time-reversal symmetry fully broken by 
B ± ), Sg(B x ,B { \) = (g(B x ,B n )) - (g(B ± » <t> Q /A,Bu)). 
Figure 2(a) shows 6g(B±, B\\) as a function of B± at 
several values of Bu , along with fits of RMT || in which 
parameters A so , t v and n± have been set by a single fit to 
the -B|| = data. The low-field dependence of 5g(Q, Bn) 
on B\\ (Fig. 2(b)) then allows the remaining parameter, 
v so , to be estimated as described below. 

Besides ez (which is calculated using g = —0.44 
rather than fit), parallel field combined with SO cou- 
pling introduces an additional new energy scale, = 

J2i j=i 2 jr~3T~' wnere «z is a dot-dependent con- 
stant and Z12 are the components of a unit vector along 
Bn H. Because orbital effects of Bii on Sg(B±, Bn) dom- 
inate at large Bn , must instead be estimated from 
RMT fits of var(g) with already-broken time reversal 



symmetry, which is unaffected by orbital coupling 24 . 

The RMT formulation [pfl i s invariant under v ao — ► 
r/v SO: where r = L1/L2 ]25|], and gives an extremal 
value of Sg(0, B^) at v so — yr. As a consequence, fits to 
Sg(0, Bn) cannot distinguish between v so and rjv so . As 



shown in Fig. 2(b), data for the 8/xm 2 dot (r ~ 2) are 
consistent with 1 < v so < 2 and appear best fit to the 
extremal value, v so ~ 1.4. Values of v so that differ from 
one indicate that both Rashba and Dresselhaus terms are 
significant, which is consistent with 2D data taken on the 
same material fiof . 

Using v so = 1.4 and values of A so , t v , and k± from the 
B|| = fit, RMT predictions for Sg(B±, Bn) agree well 
with experiment up to about B\\ ~ 0.2 T (Fig. 2(a)), 
showing a crossover from AL to WL. For higher parallel 
fields, however, experimental r5g's are suppressed relative 
to RMT predictions. By B\\ ~ 2T, WL has vanished 
in all dots (Fig. 2(c)) while RMT predicts significant re- 
maining WL at large Bn. The full range of Sg(0, Bn) for 
the three dots is shown in Fig. 2(c). The center-gated 
(5.6 /im 2 ) dot and the small (1.2 /im 2 ) dot show WL for 
all Bn , and a similar suppression of WL above B\\ ~ 2 T. 

One would expect WL/AL to vanish once orbital ef- 
fects of Bu break time reversal symmetry. Following 
Ref. (FJ), we account for this with a suppression fac- 
tor / F j(B||) = (1+r-J/r-i)- 1 , where t~\ ~ aBf + bB^ 
and assume that the combined effects of SO coupling 
and flux threading by Bn can be written as a product, 
Sg{0,B l{ ) = r5 5flA/T (0,B,|) • f F j{B\\). The B\ term re- 
fleets surface roughness or dopant inhomogeneities; the 
B^ term reflects the asymmetry of the quantum well. We 
consider fits taking a as a fit parameter (ai, Table I) with 
b = 1.410 8 s -1 T -6 fixed, obtained from self-consistent 
simulations |p6| , or allowing both a and b to be fit pa- 
rameters (02 and 62, Table I). Figure 2(c) shows that 
allowing both to be free is only significant for the (un- 
usually shaped) center-gated dot; for the small and large 
dots, the single-parameter (a) fit gives good quantitative 
agreement. 

We next consider the effects of temperature and de- 
phasing. We find that increased temperature reduces the 
overall magnitude of Sg and also suppresses AL com- 
pared to WL, causing AL at 300 mK to become WL 
by 1.5 K (maximum of Sg(B±,0) at B± = becomes 
minimum) in the 8 /im 2 dot (Fig. 3a). Fits of RMT 
to Sg(B±,0) yield A so values that are roughly indepen- 
dent of temperature (Fig. 3b), consistent with 2D results 
and T v values that decrease with increasing temper- 
ature. Dephasing is well described by the empirical form 
(^[ns])" 1 ~ 7.5T[K] + 2.5(T[K]) 2 , consistent with previ- 
ous measurements in low-SO dots p7|] . As temperature 
increases, long trajectories that allow large amounts of 
spin rotations are being cut off by the decreasing t v and 
the AL peak is diminished, as observed. 

Finally, we demonstrate in situ control of the SO cou- 
pling using a center-gated dot. Figure 4 shows the ob- 
served crossover from AL to WL as the gate voltage V g 
is tuned from +0.2 V to -IV. At V g = -IV, electrons 
beneath the center gate are fully depleted producing a 
dot of area 5.8 /itm 2 which shows WL. In the range of 
V g > —0.3 V, the region under the gate is not fully de- 
pleted and the amount of AL is controlled by modifying 
the density under the gate. Note that for Vg > 0V the 
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FIG. 4: Difference of average conductance (g) from its value at 
B± = as a function of B± for various center gate voltages V g in 
the center-gated dot (squares), along with fits to RMT fa]. Good 
fits are obtained though the theory assumes homogeneous SO cou- 
pling. Error bars are the size of the squares. Inset: A so and km as 
a function of V g extracted from RMT fits, see text. 

AL peak is larger than in the ungated 8/xm 2 dot. We 
interpret this enhancement not as a removal of the SO 
suppression due to an inhomogeneous SO coupling |p8| , 
which would enhance AL in dots with L/X so -C 1 (not 



the case for the 8/xm 2 dot), but rather as the result of 
increased SO coupling in the higher-density region under 
the gate when V g > V. 

One may wish to use the evolution of WL/AL as a func- 
tion of V g to extract SO parameters for the region under 
the gate. To do so, the dependence may be ascribed to 
either a gate-dependent A so or to a gate-dependence of 
a new parameter «|| = e|°/(((Li/Ai) 2 + (L2/A 2 ) 2 )el°). 
Both options give equally good agreement with the data 
(fits in Fig. 4 assume X so (V g )), including the parallel field 
dependence (not shown). Resulting values for A so or 
Ky (assuming the other fixed) are shown in the inset in 
Fig. 4. We note that the 2D samples from the same wafer 
did not show gate- voltage dependent SO parameters jio) . 
However, in the 2D case a cubic Dresselhaus term that 
is not included in the RMT of Ref. || was significant. 
For this reason, fits using H might show \ so {Vg) though 
the 2D case did not. Further investigation of the gate 
dependence of SO coupling in dots will be the subject of 
future work. 
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